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We report appearance of the net magnetization in Mn-based antiperovskite compounds as a result 
of the external strain gradient {flexomagnetic effect). In particular, we describe the mechanism of 
the magnetization induction in the MnsGaN at the atomic level in terms of the behavior of the 
local magnetic moments (LMM) of the Mn atoms. We show that the flexomagnetic effect is linear 
and results from the non-uniformity of the strain, i.e. it is absent not only in the ground state but 
also when the applied external strain is uniform. We estimate the flexomagnetic coefficient to be 
1.95 ^bA . We show that at the moderate values of the strain gradient (~ 0.1%) the flexomagnetic 
contribution is the only non-vanishing input to the induced magnetization. 

PACS numbers: 



Magneto-mechanical coupling in crystals has many 
practical applications, such as in sensors, magnetic 
recording devices, etc. The correlation between external 
strain and induced magnetization is in principal signif- 
icant for the systems with different dimensionality, i.e. 
bulk, thin film, and nano configurations. Yet, the cor- 
relation between the strain gradient and induced mag- 
netization is especially important in nanostructures and 
thin film heterostructures because of the large surface- 
to-volume ratio which may result in a large surface ten- 
sion due to the structural distortions caused by lattice 
mismatch, external stress, etc. For example the strain 
gradient may play a significant role when thin film is 
epitaxially grown on a substrate with slightly different 
lattice parameters. On the other hand, in problems con- 
cerning bulk structures the strain gradient is small and 
as a result has only negligible contribution. This is so 
because of the dimensional scaling inherent in the very 
definition of the strain gradient, i.e. the decrease of the 
characteristic length of the system, x results in increase 
of the Ax/x ratio. Fig. [1] schematically shows possible 
geometries of the systems under external strain gradient: 
bent nano- wire and nano-pill grown on a substrate. The 
"in-plane" and "out-of-plane" orientations of the induced 
magnetization shown on Fig. [1] will be explained later in 
the text. 

The magneto-mechanical coupling is phenomenologi- 
cally described by adding additional terms to the ex- 
pression for the thermodynamic potential, i.e. free en- 
ergy, proportional to the product of the magnetic field 
component and the conjugate terms involving mechani- 
cal strain: 
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where Xijk is the piezomagnetic tensor, Hi is the i- 
th component of the magnetic field (i = x,y,z), o~jk 
is the elastic stress tensor, /lijk is the magneto-elastic 
tensor, is the 4-rank tensor [flexomagnetic tensor). 
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is the strain gradient. For example, piezomag- 



netic (magnetostrictive) properties of certain antiferro- 
magnetic (AFM) materials are reflected in a term linear 
both in the magnetic field and in the elastic stress tensor 
For the bulk structures the last term on the r.h.s. of 
(PQ) is usually omitted because of its negligible contribu- 
tion, yet in nanostructures and/or thin film heterostruc- 
tures it plays a significant role. 

By taking partial derivative of (pp) w.r.t. magnetic field 
component, Hi we get the net magnetization in the sys- 
tem, which may be linear w.r.t. the strain (i.e. piezo- 
magnetic effect), quadratic (second order magneto-elastic 
effect), and proportional to the strain gradient: 



Mi = XijkO~jk + Vijk&jk + v ij 
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It is important to emphasize that the symmetries of 
the Xijk and fiijk tensors are different from the sym- 
metry of the Uijki tensor in ([2]), which means that the 
induced magnetization will have three contributions in 
principal distinguishable by symmetry arguments only. 
Moreover, it is possible that because of the crystal and 
magnetic symmetry the linear piezomagnetism will be 
absent in the system with the nonzero strain gradient 
induced magnetization. We define the contribution to 
the Mi from the third term on the r.h.s. of Eq. ([2]) as 
flexomagnetic (FIM) effect, i.e. strain gradient induced 
magnetization. 

To the best of our knowledge the term flexomagneto- 
electric effect was first coined by Bobylev and Pikin in 
their study of the correlation between elastic and electro- 
magnetic properties of nematic liquid crystals [2|. As a 
sidenote we emphasize here that in their work Bobylev 
and Pikin discussed the "reverse" flexomagnetoelectric- 
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ity, i.e. the re-orientation (which they called flex) of 
the molecules under external electric and magnetic fields. 
Much work was done in the past to investigate the cor- 
relation between electric polarization and external strain 
gradient, at both experimental and theoretical level Q' 
U' @|. Yet, the correlation between magnetic behavior 
of the system and the gradient of the strain was not the 
subject of mainstream research. One of the reasons for 
the insufficient study on this matter is the obvious com- 
plexity of the problem, - as opposed to the flexoelectric 
effect, where the electric polarization directly correlates 
with the atomic displacements, the F1M effect is indirect, 
i.e. it results from the re-orientation of the atomic spins 
following the atomic displacements (because of the ex- 
change interaction). As a result, for the F1M effect one 
has to consider not only the crystal but also the mag- 
netic structure and symmetry. Systems of interest must 
satisfy certain conditions, such as they have to be non- 
magnetic in ground state, and at the same time they 
have to exhibit strong magneto-elastic coupling. Com- 
paratively well known manifestation of these properties 
is the piezomagnetic effect, i.e. induction of a sponta- 
neous magnetic moment in the system under mechani- 
cal strain. In what follows we present the results of our 
study of the external strain gradient induced magnetiza- 
tion mechanism. 

To understand the mechanism of flexomagnetism at 
the atomic level we perform first principles study for 
the MnaGaN under strain gradient. The choice of the 
material is based on the intriguing magneto-mechanical 
coupling mechanism in MnaGaN. In its ground state 
MnaGaN is non-magnetic with non-collinear T 5g struc- 
ture (in the classification of Bertaut et al. [6]), i.e. the 
Mn LMMs on the (111) plane form clockwise or coun- 
terclockwise configurations, such that the spin moments 
in the plane are compensating each other. The atoms of 
Mn, Ga and N form an antiperovskite (AP) crystal struc- 
ture, and the lattice constant of the primitive 5-atom cell 
of Mn 3 GaN is 3.86 A. Fig. [2] shows the unit cell of the 
MnaGaN in the ground state. 

In our recent work [7] we have shown that the appli- 
cation of the external stress to the Mn-based AP com- 
pounds results in appearance of the non-zero magnetic 
moment, i.e. these compounds exhibit piezomagnetic 
properties. This can be explained from the symmetry 
viewpoint, i.e. as a result of the applied biaxial strain 
the symmetry of the system reduces from the trigonal 
space group PSlm to the orthorhombic Vm'm'm ferro- 
magnetic space group. Because of this transition some 
of the symmetry operations are not compatible anymore 
with the new structure. The appearance of the net mag- 
netization in the system under external stress is due to 
the rotation of the LMMs of the Mn atoms from their 
equilibrium positions. After rotation the LMMs in the 
(lll)-plane become inequivalent and do not compensate 
each other anymore. The piezomagnetic effect is linear 
and magnetization reversal is potentially possible upon 
reversal of the sign of the strain (compressive to tensile or 




FIG. 1: Schematic picture of the systems under strain gradi- 
ent with possible orientations of the induced magnetization. 



vice versa). In the present work we examine the magnetic 
behavior of the M^GaN under strain gradient (flexure), 
i.e. we look at the third term on the r.h.s. of Eq. ([2|). 




FIG. 2: (color online) Ground state of the antiperovskite 
MnsGaN unit cell: non-collinear T 59 magnetic structure. Lo- 
cal magnetic moments are shown by blue arrown on Mn 
atoms. 

We use projector augmented wave (PAW) method by 
Blochl |3L implementation of PAW by G. Kresse and D. 
Joubert [9[ in the Vienna ab initio simulation package 
(VASP) code within a Perdew-Burke-Ernzerhof (PBE) 
generalized gradient approximation [10[ of the density 
functional theory (DFT). We use a 3x12x6 k-point sam- 
pling and the Blochl's tetrahedron integration method 
jllj . We set the plane- wave cut-off energy to 300 eV and 
we choose the convergence criteria for energy of 10 ~ 5 eV. 

Accurate electronic structure calculations require pe- 
riodic boundary condition. The strain gradient breaks 
the periodicity of the primitive unit cell, therefore to 
construct a model with translational symmetry we have 
to take larger cell. For MnaGaN the smallest possible 
configuration to retain the translational symmetry un- 
der external strain gradient consists of the 8 primitive 
cells of MnaGaN (4x2x1 cell configuration). To simulate 
the strain gradient we introduce small relative atomic 
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Step 


Mn (A) 


Ga (A) 


Flex (%) 


1 


0.01868 


0.03736 


0.242 


2 


0.03736 


0.07442 


0.484 


3 


0.05589 


0.11179 


0.724 


4 


0.07450 


0.14900 


0.965 


5 


0.09264 


0.18528 


1.200 


6 


0.11132 


0.22234 


1.442 



TABLE I: Atomic shifts of Mn and Ga atoms and strain gra- 
dients: Mn24GasNs. 



shifts of Ga and Mn atoms in a way, which splits our 
40-atom cell in 4 domains (see Fig. [3j top panel). In 
each of the four domains we have strain gradient, Aa/a 
(schematically shown on the bottom panel of the Fig. [3]), 
and 4-domain configuration satisfies the periodic bound- 
ary condition. The initial orientation of the LMMs forms 
non-collinear T 5g structure. In our model we have "flat" 
triangular LMM configuration (Fig. [3l bottom panel), 
while on Fig. [2] LMMs are on the (111) plane. Both mag- 
netic orientations are energetically indistinguishable, and 
we pick the " flat" system for simplicity. Table [J sum- 
marizes the values for the atomic shifts and the strain 
gradients for the Mn24GagN8 cell. 
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FIG. 3: (color online) Top panel: Mn24GasNg cell: 4-domain 
structure. Mn atoms - blue, Ga atoms - grey, N atoms - dark 
red (almost invisible behind Mn atoms). The directions of 
the atomic shifts are schematically shown by black arrows. 
Bottom panel: Single domain in Mn24GasNs under strain 
gradient. 



Next we relax the LMMs of the Mn atoms but we 
keep the atomic positions fixed to make sure that the 
crystal structure does not relax back to the unstrained 



ground state. Fig. [H shows results of our calculations 
for the magnetization per Mn atom as a function of the 
strain gradient. There are two different contributions: 
the blue line with circles represents out-of-plane mag- 
netization (multiplied by 10) which is a direct result of 
the flexure. If the strain applied to the cell is uniform 
then from the symmetry arguments it is clear that only 
in-plane magnetization will appear. The black line with 
squares represents in-plane non-linear contribution to the 
magnetization which results from the in-plane rotations 
of the LMMs of the Mn atoms (shown on the bottom 
panel of the Fig. [3] by black arrows). Important fea- 
ture of these two mechanisms is that at the moderate 
values of the strain gradient the dominant contribution 
comes from the linear flexomagnetic magnetization, while 
the non-linear in-plane contribution is vanishingly small. 
This correlation changes at the higher values of the gradi- 
ent but the dominant nature of the linear contribution at 
the moderate gradient values is important for the prac- 
tical applications where the regular values of the strain 
gradient are of the order of ~ 0.1 0.2%. 
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FIG. 4: (color online) Linear out-of-plane (xlO) and non- 
linear in-plane induced net magnetization of Mn atom as a 
function of the strain gradient. 



To understand better the nature of these two contri- 
butions we examine the behavior of in-plane and out-of- 
plane components of the LMMs of Mn atoms within one 
domain. Figures [5] and [6] show our results. The out-of- 
plain contribution comes from the Mn atoms 3, 4 and 7 
(see Fig.[3j bottom panel). At the same time the in-plane 
rotations of the Mn LMMs (atoms 3 and 7) show distinct 
non-linear feature, which contributes to the non-linear in- 
plane magnetization shown on Fig. [H The appearance of 
the out-of-plane component is not observed if the applied 
strain is uniform. Therefore, the appearance of the lin- 
ear out-of-plane magnetization is purely strain gradient 
related effect and its mechanism is different from the one 
responsible for the in-plane magnetization induction. Be- 
low we present our phenomenological arguments on the 
mechanism of the F1M effect at the atomic level and es- 
timate the value of the F1M coefficient. 

Qualitatively, the out-of-plane rotations of the LMMs 
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FIG. 5: Out-of-plane magnetic moment vs. strain gradient. 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 



N 

CD 
C 
03 
Q_ 





Mn:2,4,8 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 






Mn: 3, 7/ 


-1.18 
-1.19 
-1.20 
-1.21 
-1.22 


Mn: 1,3,5jJ, 9 


i i i i i i i 


i 



Strain gradient (%) 

FIG. 6: In-plane magnetic moment vs. strain gradient. 



can be explained by the frustration of the spin orien- 
tations. Because of the inequivalent strain on differ- 
ent lattice sites the ground state nearest neighbor dis- 
tances between Mn atoms change in non-uniform way. 
This is schematically shown on the bottom panel of the 
Fig. [3] where for example the Mnl^Mn4 distance be- 
comes larger than the Mn4^Mn5 distance. This is a 



pure strain gradient related effect, this non-uniformity is 
not present in the ground state and in the case of the 
uniform strain. This results in the increase of the ex- 
change integral, Jfr^s) in the Heisenberg Hamiltonian 
between Mn4 and Mn5 atoms (see Fig. G2 bottom panel) 
due to the increasing overlap of their wave functions and 
at the same time in the decrease of the J(ri^) between 
Mnl and Mn4 atoms. This mechanism is responsible for 
the out-of-plane rotation of the Mn4 (and other "inside 
of the domain" Mn atoms, such as Mn3 and Mn7 on the 
Fig. bottom panel) LMM. 

We estimate the flexomagnetic coefficient at the strain 
gradient value of ~ 0.4% from Eq.© and the magnetiza- 
tion vs. strain gradient data shown on Fig. [4] as follows. 

Mfi ex — U^] 

da _ d [ Afl i_ d r ao+c-gn 1 0.004 
dx dxt a 1 dxt a -J a 3.9-10- 10 m' 

M fiex 0.002/iB , lt -I Qr,, X 

V = -gp= .004/3.9.10-i»m V ~ L95 M£ A ' 

Here ao = 3.9 • 10 _10 m is the ground state lattice con- 
stant of the MnaGaN, Mfi ex = 0.002/i# is the induced 
out-of-plane magnetization value at ~ 0.4% of the strain 
gradient. Since the magnetization is linear, the v will 
have the same value over the considered range of the 
strain gradient. 

In summary, in this paper we discussed the appear- 
ance of the net magnetization in Mn-based antiper- 
ovskite compounds (in particular, MnaGaN) under ex- 
ternal strain gradient [flexomagnetic effect). The mag- 
netization dependance on the flexure is linear. The es- 
timated flexomagnetic coefficient is 1.95 /i#A. Besides 
being of substantial theoretical value, the flexomagnetic 
effect can have interesting practical applications, such 
as in electrical control of magnetization in memory cells 
if the structural distortion in the flexomagnetic phase 
is induced and can be controlled by the external elec- 
tric field (i.e. by forming heterostructure with ferroelec- 
tric/piezoelectric compounds). 
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